Kaon Interferometry: A Sensitive Probe of the QCD Equation of State? 
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We calculate the kaon interferometry radius parameters for 
high energy heavy ion collisions, assuming a first order phase 
transition from a thermalized Quark-Gluon-Plasma (QGP) to 
a gas of hadrons. At high transverse momenta Kt ~ 1 GeV/c 
direct emission from the phase boundary becomes important, 
the emission duration signal, i.e., the -Rout/flWe ratio, and 
its sensitivity to T c (and thus to the latent heat) are en- 
larged. The QGP+hadronic rescattering transport model cal- 
culations do not yield unusually large radii (Ri < 9fm). Fi- 
nite momentum resolution effects have a strong impact on the 
extracted interferometry parameters (Ri and A) as well as on 
the ratio i? ou t/-Rsidc- 

In this Letter we present predictions for the kaon in- 
terferometry measurements in Au+Au collisions at the 
Relativistic Heavy Ion Collider (RHIC) (nucleon-nucleon 
center-of-mass energies up to y/s NN — 200 GeV). Cor- 
relations of identical kaon pairs as well as Bose-Einstein 
correlations in general represent an important tool for 
the understanding of the space-time dynamics in multi- 
particle production processes In the particular case 
of relativistic heavy ion collisions, one important goal is 
to prove the existence of a phase transition from quark- 
gluon matter to hadrons as predicted by QCD lattice 
calculations for high temperatures. Moreover, the prop- 
erties of that phase transition, for example the critical 
temperature T c or its latent heat, are of great interest. 
For a first order phase transition, the associated large 
hadronization time was predicted to lead to unusually 
large interferometry radii ||-fl| . The mixed phase of co- 
existing hadrons and partons prolongs the emission dura- 
tion of particles from the phase boundary because of the 
large latent heat. This should then lead to an increase of 
the effective source size, in particular in the outward di- 
rection, i.e., parallel to the transverse pair velocity. This 
phenomenon was also expected to depend on the critical 
temperature T c , the latent heat or the specific entropy 
of the quark-gluon phase (for recent reviews, see, e.g., 
^,[1). However, as demonstrated recently 0, the sub- 
sequent hadronic rescattering phase following hadroniza- 
tion from a thermalized QGP does not only modify the 
pion interferometry radii but even dominates them. Late 
numerous soft collisions in the hadronic phase diminish 
and also alter qualitatively the particular dependencies 
on the QGP-properties 

Here, we investigate what can be learned from the in- 
terferometry analysis of kaons. The motivation to extend 
the interferometry analysis from pions to kaons (which 
are measured soon at RHIC) is provided by several as- 



pects, (i) Kaons are expected to be less contaminated 
by resonance decays compared to pions (ii) In the 

case of neutral kaon correlations, two-particle Coulomb 
interactions do not distort the Bose-Einstein correlations 
(while non-negligible strong final state KK interactions 
persist due to the near threshold resonances a0(980) 
and f0(980) jnj). (iii) The kaon density is consider- 
ably smaller than the pion density. The pion multiplic- 
ity itself has increased by approximately 70% from SPS 
(^iv = 17.4 GeV) to RHIC (^ NN = 130 GeV) 0] 
whereas the interferometry radii are almost the same 
fli^pm . Hence, higher multiparticle correlation effects 
that might play a role for pions, should be of minor im- 
portance for kaons. The comparison of kaon and pion 
correlation functions will provide a test of the presently 
applied two-particle correlation formalism. (iv) The 
strangeness distillation mechanism fl4| might further in- 
crease the time delay signature -R ut/-Rsidc- Kaon evap- 
oration could lead to strong temporal emission asymme- 
tries between kaons and antikaons [Hj] , thus probing the 
latent heat of the phase transition. 

We will show in the following that for kaons, as for 
the pions, the bulk properties of the two-particle corre- 
lation functions are dominated by a long-lived hadronic 
rescattering phase. Thus, the interferometry radii appear 
to depend only weakly on the precise properties of the 
QGP, such as the thermalization time Ti, T c , the latent 
heat or the specific entropy of the QGP. However, we will 
demonstrate that this sensitivity is considerably enlarged 
at high transverse momenta Kt ~ 1 GeV/c. In this kine- 
matic regime, direct emission from the phase boundary 
becomes important (~ 30%) allowing us to inspect the 
prehadronic phase with less distortions. Finally, we will 
calculate explicitly the correlation functions and extract 
the corresponding parameters with and without taking fi- 
nite momentum resolution effects into account. A strong 
impact on the radii, the A intercept parameter and the 
-Rout/-Rsidc ratio are observed; they all decrease, depend- 
ing on Kt, if a finite momentum resolution is taken into 
account. 

The calculations are performed within the framework 
of a relativistic transport model that describes the initial 
dense phase of a QGP by means of ideal hydrodynamics 
employing a bag model equation of state that exhibits a 
first order phase transition j§,[l6]]. Hence, we focus on a 
phase transition in local equilibrium, proceeding through 
the formation of a mixed phase. Smaller radii and emis- 
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sion times may result for a crossover jl 17 or for a rapid 
out-of-equilibrium phase transition similar to spinodal 
decomposition |lq| . Cylindrically symmetric transverse 
expansion and longitudinally boost-invariant scaling flow 
are assumed [^| Jl^Jl^| . This approximation should be rea- 
sonable for central collisions at high energy, and around 
midrapidity. The model reproduces the measured pt- 
spectra and rapidity densities of a variety of hadrons at 
V$nn = 17.4 GeV (CERN-SPS energy), when assuming 
the standard thermalization (proper) time Tj = 1 fm/c, 
and an entropy per net baryon ratio of s/ ps = 45 Jl9| , p0| . 
Due to the higher density at midrapidity, thermalization 
may be faster at BNL-RHIC energies - here we assume 
t, = 0.6 fm/c and s/pb = 200. (With these initial con- 
ditions results on the multiplicity, the transverse energy, 
the p<r-distribution of charged hadrons, and the p/p ra- 
tio are described quite well [ ^9[p0| ].) The later hadronic 
phase is modeled via microscopic transport that allows 
us to calculate the so-called freeze-out, i.e., the time and 
coordinate space points of the last strong interactions 
of an individual particle species, rather than applying 
a freeze-out prescription as necessary in the pure hy- 
drodynamic approach. Here, we employ a semi-classical 
transport model that treats each particular hadronic re- 
action channel (formation and decay of hadronic reso- 
nance states and 2 — > n scattering) explicitly pl| . The 
transition at hadronization is performed by matching the 
energy-momentum tensors and conserved currents of the 
hydrodynamic solution and of the microscopic transport 
model, respectively (for details, see Q). The micro- 
scopic model propagates each individual hadron along a 
classical trajectory, and performs 2 — ► n and 1 — * m pro- 
cesses stochastically. Meson-meson and meson-baryon 
cross sections are modeled via resonance excitation and 
also contain an elastic contribution. All resonance prop- 
erties are taken from fz^| . The ttK cross section for ex- 
ample is either elastic or is dominated by the if* (892), 
with additional contributions from higher energy states. 
In this way, a good description of elastic and total kaon 
cross sections in vacuum is obtained . Medium effects 
on the hadron properties, as for example recently studied 
by hydrodynamical calculations employing a chiral equa- 
tion of state |L7]], are presently neglected. For further 
details of this dynamical two-phase transport model, we 
refer to refs. 

For the correlation analysis, a coordinate system is 
used in which the long axis (z) is parallel to the beam 
axis, where the out direction is parallel to the trans- 
verse momentum vector Kt = (pit + P2t)/2 of the 
pair, and the side direction is perpendicular to both. 
Thus, -Rout probes the spatial and temporal extension 
of the source while -R s idc only probes the spatial exten- 
sion. Thus the ratio -R ou t /-Rsidc gives a measure of the 
emission duration (see also eqs.(l)-(3) and discussion be- 
low). It has been suggested that the ratio -R ut/-R s ide 
should increase strongly once the initial entropy density 
Si becomes substantially larger than that of the hadronic 



gas at T c g. The Gaussian radius parameters are ob- 
tained from a saddle-point integration over the classi- 
cal phase space distribution of the hadrons at freeze-out 
that is identified with the Wigner density of the source, 
S(x,K) @!1§. 

i? b 2 idc (K T )Hy 2 }(K T ), 



R? 



„ut(K T ) 

long( K T) 



{{x - p t t) 2 )(K T ) = (x 2 + fit 2 - 2(3 t xt) 
((5-Ai) 2 )(K T ), 

with ^(Kt) — x^ — (x a ')(Kt) being the space- 
time coordinates relative to the momentum depen- 
dent effective source centers. The average in (|l|)-(||) 
is taken over the emission function, i.e. (/) (K) = 
/ d 4 xf(x)S(x,K)/ J d 4 xS(x,K) with K = (E K ,K). In 
the osl system j3 = , 0, /?;), where (3 = "K/Ek and 
Ek = V m 2 + K 2 . Below, we cut on midrapidity kaons 
(|y| < 0.5). In the absence of x-t correlations, i.e. in par- 
ticular at small Kt, a large duration of emission At = 
\J{t 2 ) increases -R ou t relative to i? s ide [&fl|- F° r strong 
(positive) x-i-correlations or large spatial anisotropics 
in out- and side-direction ((y 2 ) > (x 2 )), in principle 

-Rout < -Rsidc COuld follow. 

The absolute values of the kaon radii determined by the 
above expressions (|l])-(||) are considerably smaller than 
the pion radii, especially at low Kt- These pion radii 
are larger by a factor of two at low Kt (< 400 MeV) 
while at higher Kt the values become similar. This 
is due to the resonance source character of pions. Mi- 
croscopic transport calculations show that at SPS en- 
ergies (v^ATAr — 17.4 GeV) about 80% of the pions are 
emitted from various resonances J25|. This leads to a 
strong substructure of the freeze-out distributions p5[ , 
e.g. strongly non-Gaussian tails. The ratio -R ut/i?sidc 
for kaons is shown in Fig. 1. Due to identical freeze-out 
distributions the correlation parameters are the same for 
kaons and antikaons. The bag parameter B is varied 
from 380 MeV/fm 3 to 720 MeV/fm 3 , (i.e., the latent 
heat changes by ~ 4-B), corresponding to critical tem- 
peratures of T c ~ 160 MeV and T c ~ 200 MeV, respec- 
tively. A change of T c implies a variation of the lon- 
gitudinal and transverse flow profiles on the hadroniza- 
tion hypersurface (which is the initial condition for the 
subsequent hadronic rescattering stage). i? ut /-Rsidc is 
smaller for kaons than for pions at the same small Kt 
because of their larger mass (leading to smaller veloci- 
ties which reduce the temporal contribution to -R ou t in 
eq. (|) ((3f(t 2 )) |§). The ratios reach a value of 1.5 at 
Kt ~ and Kt ~ rriK for the rather rapidly rising 
pions 0| and the gradually growing kaons, respectively. 
Model calculations solely based on hadronic degrees of 
freedom generally predict smaller ratios p7|j9],p8t. The 
sensitivity of the value of -R ou t /-Rside to the critical tem- 
perature T c increases strongly with Kt- Higher T c speeds 
up hadronization but on the other hand prolongs the dis- 
sipative hadronic phase that dominates the radii. More- 
over, in the lower T c case, direct emission and immediate 
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freeze-out from the phase boundary becomes important 
at large Kt (~ lGeV/c). High- Kt kaons are strongly 
correlated with early mean emission times. The reso- 
nance contribution for the kaons is still quite large, de- 
creasing with K T from 70 to 50% for T c ~ 160 MeV. How- 
ever, most of these kaons are from if* (892) decays with 
the K* having a moderate lifetime of r w 4 fm/c. Elastic 
scatterings prior to freeze-out contribute on the order of 
20%. The direct emission from the phase boundary, i.e., 
the kaon did not suffer further collisions in the hadron 
gas after the particle had hadronized, increases strongly 
(approximately linearly with Kt) for T c ~ 160 MeV up 
to 30% at K T = 1 GeV/c. For the higher T c (~ 200 MeV) 
hadronization is earlier. The hadronic phase lasts longer 
and the system is rather opaque for direct emission. This 
direct emission component is not present in pure ideal 
hydrodynamical calculations (e.g. ]26|]) for which all par- 
ticles, also at high Kt, are in (local) thermodynamical 
equilibrium. Thus, there is no possibility for direct emis- 
sion from the phase boundary and escaping the hadronic 
phase unperturbed. 
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FIG. 1. -Rout/-Rsid c as obtained from eqs. (1) and (2) for 
kaons at RHIC (^/s NN = 200 GeV) (full symbols) and at SPS 
(y/s NN = 17.4 GeV) (open symbols), as a function of Kt 
for critical temperatures T c ~ 160 MeV and T c ~ 200 MeV, 
respectively. The crosses show pions for comparison. The 
lines are to guide the eye. 

Finally, we calculate the correlation parameters (shown 
in Fig. 2) by performing a \ 2 fit of the three-dimensional 
correlation function C2((fout, <7side, QWg) to a Gaussian as 

C 2 (g , 9s, ffl) = 1 + Aexp(-«£i£ - q 2 s R 2 s - q 2 R 2 ) . (4) 

The correlation functions are calculated from the phase 
space distributions of kaons at freeze-out using the cor- 
relation after burner by Pratt p|,p3| . It is assumed that 
the particles are emitted from the large system indepen- 
dently, which allows to factorize the iV-boson production 
amplitude into N one-boson amplitudes A(x). Then, the 



emission function is computed as the Wigner transform 
S(x,K) = Jd 4 ye i y K A*(x + y/2)A(x-y/2). The two- 
boson correlation function is given by 

C 2 (p,q)-1 = 

/ d 4 xS(x, K) / d 4 yS(y, K) exp (2ik ■ (x - y)) 



fd*xS(x,p)fd*xS(y t <0 
J d 4 xS(x, K) / d 4 yS(y, K) exp (2ik ■ (x 



y)) 



\ J d 4 xS(x,K)\ 2 



(5) 



E p — Eg. The 



where 2K = p + q, 2k = p q, and 2k° 
second line in (^) holds in the limit where the width of 
the correlation function is small such that p ~ q ~ K 
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FIG. 2. Kaon correlation parameters as obtained from a % 2 
fit of C'2 (eq. (5)) to the Gaussian ansatz (eq. (4)) for Au+Au 
collisions at RHIC (v^jvjv ~ 200 GeV). The top and bottom 
panels show the results for T c ~ 160 MeV and T c ~ 200 MeV, 
respectively. R ou t (circles) and -R s id c (squares) are shown on 
the left while -Ri on g (diamonds) and A TO (triangles) are shown 
on the right. Full and open symbols correspond to calcula- 
tions with and without taking momentum resolution effects 
into account, respectively. 

Given a model for a chaotic source described by 
5(x, K), such as the transport model described above, 
eq. (|^) can be employed to compute the correlation func- 
tion. The expressions based on an Gaussian 
ansatz, yield larger values for the pion radii than per- 
forming a fit to the correlation functions. For the kaon 
transverse radii, similar results are obtained with both 
methods. Only i?i ong is larger if determined by (||) as 
compared to the fitting result. This is due to the non- 
Gaussian contribution arising from the longitudinal ex- 
pansion dynamics that is similar for pions and kaons |p8| . 
Kaons are better candidates for the Gaussian expressions. 
There are fewer resonance decays into kaons compared to 
pions, and long-lived resonances do not play as important 
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a role. 

For T c ~ 200 MeV, i? ut is only approximately 1 fm 
larger than in the T c ~ 160 MeV case. This reflects a 
fact already known from the pions. Higher T c leads to an 
earlier hadronization, thus causing a prolonged hadronic 
phase. When taking finite momentum resolution (f.m.r.) 
into account, the true particle momentum p obtains an 
additional random component. This random component 
is assumed to be Gaussian with a width Sp. The relative 
momenta of pairs are then calculated from these modi- 
fied momenta. However, the correlator is calculated with 
the true relative momentum. While i? ut remains con- 
stant or even slightly increases with Kt when calculated 
without f.m.r., it drops if a f.m.r. of w 2% of the center 
of each Kt bin is considered, a value assumed for the 
STAR detector Accordingly, the differences with 

and without f.m.r. increase with Kt- The f.m.r. leads 
to smaller radii. i? ut is strongly reduced while i? s ide 
shows a moderate reduction. Thus, the -R ou t/-R S idc ratio 
is considerably reduced through the f.m.r.. For exam- 
ple, in the T c ~ 200 MeV case, it is reduced from 1.8 to 
1.35. However, it is always larger than one. A proper 
correction for f.m.r. (STAR and other experiments ]29| ] 
do correct) is needed to maintain the sensitivity to the 
QGP-properties. The A parameter is roughly constant as 
function of Kt for Sp/p — but it decreases rapidly with 
a, f.m.r.. The correlation strength is transported to larger 
q values by the f.m.r. effects. 

We have calculated kaon correlation parameters for 
Au+Au collisions at RHIC energies, assuming a first- 
order phase transition from a thermalized QGP to a gas 
of hadrons. No unusually large radii are seen (Ri < 
9 fm) . A strong direct emission component from the 
phase boundary is found at high transverse momenta 
(Kt ~ lGeV/c) where also the sensitivity to the crit- 
ical temperature, the latent heat and specific entropy 
of the QGP is enlarged. Finite momentum resolution 
effects reduce the true parameters (Ri, A) and the ra- 
tio i?out/-Rsidc substantially. Kaon results from RHIC at 
high Kt will provide an excellent probe of the space-time 
dynamics close to the phase-boundary and to the prop- 
erties of this prehadronic state, possibly an equilibrated 
Quark-Gluon-Plasma. 
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